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, a fermionic Chern-Simons approach we study the description of the system's 



field at and near an even denominator Landau level filling fraction. Using 



^ \ low energy excitations within an extension of Landau's Fermi liquid theory. 

^ \ We calculate perturbatively the effective mass and the quasi-particle interac- 

Q tion function characterizing this description. We find that at an even denom- 



inator filling fraction the fermion's effective mass diverges logarithmically at 



the Fermi level, and argue that this divergence allows for an exact calculation 



of the energy gaps of the fractional quantized Hall states asymptotically ap- 
proaching these filling fractions. We find that the quasi-particle interaction 
function approaches a delta function. This singular behavior leads to a cance- 
lation of the diverging effective mass from the long wavelength low frequency 
linear response functions at even denominator filling fractions. 
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I. INTRODUCTION 



The Fermionic Chern-Simons theory has been extensively used in the last few years to 
describe the physics of a two dimensional electron gas in a partially filled Landau level [[I 
|J. In that theory, the problem of electrostatically interacting electrons in a strong magnetic 
field is converted, using an exact transformation, to a problem of electrostatically interacting 
"transformed fermions" which also interact with a Chern-Simons gauge field. 

Being still unsolvable, the transformed fermion problem is studied by approximations. 
The simplest approximation, mean field theory, already gives interesting physical results. 
For example, it predicts the existence of stable quantized Hall states at filling fractions of 
the form j^^, where p is an arbitrary integer, and is an even positive integer, equal to 
the number of quanta of Chern-Simons flux attached to each fermion. As was first noted by 
Jain 0, filling fractions of the above form, with = 2, are the most prominent fractional 
Hall states seen in the lowest Landau level. (Jain's "composite fermions" may be regarded as 
equivalent to the transformed fermions, projected onto the lowest electronic Landau level.) 
By contrast, at the even denominator fractions 1/0, the fermion Chern-Simons mean field 
theory predicts a Fermi-liquid state, with no energy gap and no quantized Hall effect |J [[|. 

For a description of the excitation spectrum and of other quantities of interest, it is of 
course necessary to go beyond mean field theory. Specifically, it is necessary to consider, 
using perturbation theory, the effects of fluctuations in the Chern-Simons field and of the 
two body interaction, effects which are omitted from mean field theory. Consequences of 
these perturbations were explored at some length by Halperin, Lee and Read (HLR) fl|] for 
the Fermi liquid state at v = 1/0, and also to some degree for other nearby filling fractions. 

An interesting result of the HLR analysis was the prediction that at v = 1/0 the effective 
mass m* of the transformed fermions diverges for quasi-particles at the Fermi energy. For 
the physically relevant case of Coulomb interaction between the electrons, this divergence 
was found to be quite weak, i.e., logarithmic in u, the energy relative to the Fermi energy. 
For short range electron-electron interaction HLR found a stronger divergence m* ~ cu^ 1 ^ 3 , 
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while for a longer range interaction the effective mass remains finite as 00 — > 0. 

The HLR results were based on an examination of the leading diagrams in the perturba- 
tion expansion for small <f), supplemented by heuristic arguments that the leading singularity 
might not be affected by higher orders in perturbation theory. Consequently, it is of interest 
to see whether these arguments can be made more convincing, or whether by contrast results 
contradictory to HLR may be obtained. Interest in these questions has been stimulated by 
the fact that similar effective mass divergences are found in other models in which fermions 
interact strongly with a fluctuating gauge field, including models that have been advanced 
in the context of high T c superconductivity Q- ||11|| . 

Assuming that the effective mass divergence predicted by HLR is indeed correct, one 
may ask what physically measurable quantities might reflect it. According to HLR, the 
renormalized mass should not show up in the long wavelength density and current response 
functions at zero temperature, for uj/q — > 0. This observation was confirmed, more recently, 
by analyses using a variety of theoretical techniques, including renormalization group meth- 
ods U ||, bosonization [Q and diagrammatic analyses [[J. Moreover, Kim et.al. [|], using an 
approximation that takes into account the leading divergent contribution to m*, have shown 
that there is no corresponding singularity in the density response function in the limit of 
low frequencies and long wavelengths, for any value oiuo/q. Within their approximation, the 
density response function looks just like the one obtained in the random phase approxima- 
tion of a Fermi gas with the bare mass m. The diverging effective mass may be expected to 
show up in the temperature dependence of the specific heat at low temperatures; however, 
a contribution from the low frequency density relaxation collective mode could lead to the 
same anomalous temperature dependence as the quasiparticle contribution. 

As noted by HLR, a more definitive manifestation of the diverging m* should come from 
the behavior of the energy gap at quantized Hall steps for v = ^j-p for large values of p. 
Specifically, one might expect the energy gap E g (u) to be 




(1) 
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where m*{y) is the effective mass calculated self-consistently at an energy E g {y), and AB = 
B — 2ir(f)n is the deviation of the magnetic field from its value at v = 1/0. Here we denote the 
average density of electrons by n, and we use units where h — - — 1. Until now, however, 
the conjecture ([1]) has not been checked by any explicit calculation. 

In this paper we study both the even denominator state and the way it is approached 
when the filling fraction is slightly away from it. We study the properties of the Fermi liquid 
formed at an even denominator filling fraction (e.g., 1/2), namely, the effective mass and the 
quasi-particle interaction function. We review the perturbative calculation leading to the 
divergence of the effective mass, and find that for the Coulomb case, as uj —>■ 0, the leading 
term in the effective mass is, 



<- 2 ek, 



m» = ^|lnu;| (2) 

where e is the background dielectric constant, uj is the distance from the Fermi energy and 
the Fermi wavevector hp is related to the electron density by kp = (47m) 1 / 2 . We then argue 
that Eq. (0) is the exact leading term in the limit uj — > 0. This statement is based on the use 
of Ward's identity to account for the renormalization of the vertices in the diagram leading 
to the divergence of m* . We directly calculate the energy gap at a filling factor of the form 



p 



(f>P+i 



and confirm the conjecture ([!]). Thus, for large p, the energy gap we obtain is, 

E 3 * U (3) 

9 3 /2 el H (4>p + 1) In (4>p + 1) 

where Ih = \f^§ is the magnetic length. Following our discussion of the effective mass, we 



argue that our expression for the energy gap, Eq. (Q), is exact in the p — > oo limit |L3 
In our study of the Landau quasi-particle interaction function at v = 1/0 we find that it 
includes a singular contribution that is just of the form necessary to explain the results of 
Kim et.al. 0, i.e., of the form needed to cancel the effect of a diverging effective mass on the 
zero temperature long wavelength low frequency linear response function. The expression we 
find for the Landau quasi-particle interaction function is, however, not exact. In particular, 
our approximation does not properly account for the effects of short wavelength and high 



frequency fluctuations. These fluctuations are particularly important in the limit where the 
electron's bare mass is vanishingly small, so that the cyclotron energy becomes infinite and 
the actual electron states are projected onto the lowest Landau level; our approximation 
incorrectly predicts the magnitude of the Landau interaction function in this limit. The 
short wavelength and high frequency fluctuations are not important for the effective mass 
m* and the energy gaps at p — > oo, because these quantities are determined by infra-red 
divergences in the self energy. 

Combining our analysis of the effective mass and the quasi-particle interaction function, 
we study the behavior of the chemical potential n of the composite fermions as the density 
n is varied and AB is kept fixed. We find that the chemical potential jumps discontinuously 
whenever the density corresponds to an integer fermionic filling factor p. The magnitude of 
the jump is the energy gap corresponding to the fractional quantum Hall state at v = t^j- 
As emphasized above, it reflects the divergence of m*. However, we also find that the 
chemical potential varies continuously during the course of a population of a fermions' 
Landau level. This variation is such that the overall slope of fi(n) is independent of the 
divergent contribution to m* (See Fig. (|5|)). 

In our perturbative study of the Fermion's effective mass, energy gaps and Landau in- 
teraction function, we use the RPA expression for the gauge field propagator. This use is 
then justified by our results. The singularities we study result from the long wavelength low 
frequency behavior of the gauge propagator. The latter are determined by the long wave- 
length low frequency limit of the fermion's polarization functions (the fermions irreducible 
polarization bubble). As we show, and as was shown before by several authors @ @ @, 
in that limit, the fermions' irreducible polarization functions are not renormalized by gauge 
field fluctuations. 

The outline of the paper is as follows: in Section (2) we construct a general framework 
for the discussion of the Fermi liquid formed by the composite fermions. In Section (3) 
we calculate perturbatively the parameters characterizing this Fermi liquid, and we discuss 
the regime of validity of the perturbative approach. In Section (4) we examine the linear 
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response function resulting from the parameters obtained in Sec. (3). In Section (5) we 
calculate the energy gaps for large p, and analyze the dependence of the chemical potential 
on the density for fixed AB. In Section (6) we comment on the effect of the coupling of the 
fermions to longitudinal gauge fluctuations. We conclude with a summary. 

Several details of the analysis of Section (3) are given in Appendices A and B. 

II. GENERAL DESCRIPTION OF THE FERMI LIQUID OF THE COMPOSITE 

FERMIONS 

In this section we define the model to be considered, and we give a general description 
of the Fermi liquid formed at even denominator filling fractions. We consider electrons in 
a two dimensional electron gas, subject to a magnetic field B. The electrons interact via a 
Coulomb interaction, V(r — r') = t ^l r i\ ■ Using a singular gauge transformation, the problem 
is mapped onto that of fermions subject to a magnetic field AB = B — 2im(f). (In our system 
of units, where h — - — 1, the flux quantum is 2tt.) In the case we discuss in this section, 
the electronic filling factor is 1/0, and then AB = and the fermionic filling factor, p, is 
infinite. In section (5) we discuss the case of finite AB and p. Generally, the fermions' 
Hamiltonian is 

H = ^ J d 2 rvEr+(r)[-iV + AA(r) - a(r)] 2 $(r) 

(4) 

+U d 2 r J d 2 r'[^ + (r)^/(r) - n]V(r - r')[^ + (r')^(r') - n] 

supplemented by the gauge condition V • a = and by the constraint V x a(r) = 
27r0[\I /+ (r)\[ / (r) — n], where n is the average electron density. In Eq. (|J) AA(r) = \ABz x r. 
The constraint may be used to write the electrostatic term as 2 ( 27 ^ 2 / ( ^' r I d 2 r'[V x 
a(r)]V(r — r')[V x a(r')]. The fermion's interaction with the transverse gauge field is then 
of the form —J • a — T^pa 2 where 3,p are the fermions' current and density. The Fermi 
wave- vector of the fermions is kp = 
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At the mean field approximation, the Chern-Simons gauge field a is replaced by its 
average value, and thus its interaction with the fermions is neglected. In this section we 
attempt to describe the effect of these interactions on the fermions within an extension of 
Landau's Fermi liquid theory. Landau's theory describes the low energy excitations of the 
system in terms of the occupation number function n(k) ||15|| - [17]. This function signifies 



the difference between the occupation number of the state k at the excited state and the 
corresponding number at the ground state. The energy density corresponding to a given 
function n(k) is 

£ ("< k » = 1 + \l W?l Sf s ^ V)n ^ ] (5) 

where = is the kinetic energy of a quasi-particle near the Fermi level, m* is its 
effective mass, and /(k, k') is the Landau interaction function. Note that in the presence of 
a non-zero vector potential, k is the dynamical momentum, and not the canonical one. 

The parameters m*, / include the entire effect of the interaction near the Fermi level in 
the case of a neutral fluid such as 3 He in its normal state, but they conventionally leave out 
the direct, Hartree, part of the interaction in the case of electrons in a metal. The difference 
between the two cases is in the range of the interactions being short for Helium atoms and 
long for electrons in metals. As developed in Silin's extension of Landau's theory |17|| , the 



effect of long range interactions is taken into account when linear response functions are 
calculated, by making a distinction between the externally applied driving force and total 
driving force. In the problem we consider, the Chern-Simons interaction of the fermions with 
the gauge field is long ranged. Consequently, the direct Hartree part of both that interaction 
and the electrostatic one should be separated from the rest. Thus, our construction of the 
Fermi liquid picture for the composite fermions starts by defining the energy functional ([5]) 
excluding the contribution of the direct Hartree part of both interactions. Then, we use the 
equation of motion derived from that functional to define and calculate a linear response 
function of the fermions. 

The equation of motion is conveniently written as an integral equation for the function 



v{9), defined by n(k) = 5(k — k F )u(9), where 9 is the angle between k and some reference 
direction jn|- ||17|| . The function v(9) describes the deformation of the Fermi surface in the 
direction 9. The equation of motion for v(9) is 

(a; - qv* F cos 9)^(9) - ^ q k F cos9 J d9'f{9 - 0> qjW (0') = k • F q , w (6) 

where the angle 9 is the angle between the vectors k and q; the function f{9 — 9') is the 
quasi-particle interaction function f(k, k') for two wave-vectors at the Fermi surface, at 
angles 9, 9' relative to q; F is the driving force, characterized by a wave-vector q and a 
frequency u\ and the renormalized Fermi velocity is denoted by v* F . Eq. (||) is valid for 
q <^ k F and u <ti fi. 

In Landau's Fermi liquid theory the linear response function is extracted from Eq. (^|). 
This equation relates the deformation function u(9) to the driving force F. In terms of 
v(9), the fermions' density is p qiW = k F J d9h , qLtUJ (9) while the current density is J qA , = 
/rf 2 M(fc-fc F )u k z/ qiW (^), withu k = ±+J d 2 /c7(k,k')<5(e k /-/x)^. The driving force can be 
expressed in terms of derivatives of a three-vector potential A u , as F qja , = — iuA^+iqAo^^. 
Thus, by means of Eq. (|6]), the current vector J M can be related to the vector potential A u by 
J M = U 11U A U , and the linear response function can be calculated. Generally, the linear 
response function is a matrix relating the current-density three-vector J M to the three-vector 
potential A v . Since the three- vector J M is constrained by the conservation of charge and the 
three-vector A u is constrained by a gauge condition V • A = 0, the linear response function 
can in fact be described by a 2 x 2 matrix ILj^q, ou), with the two indices taking the values 
(for the time component), and 1 (for the transverse component). 

Now we turn to discuss the effect of the Hartree term. The function II relates the current 
J M to the total vector potential acting on the fermions. This vector potential is composed of 
an externally applied part and a part induced by the fermion current and density. A fermion 
current-density vector J M induces a vector potential A™ d = V^J^ where the matrix V is 
given by, 



S 



V 



( 2-ire 2 2-k%4> \ 
eq q 

V 1 ) 



(7) 



The diagonal term is the electrostatic potential induced by a charge density. The off diagonal 
terms describe the vector potential induced by the fictitious flux tubes attached to the 
fermions. 

We can now define the matrix K, relating the current-density vector to the externally 
applied vector potential A e v xt by J M = K^A^K Since A v = A ext + A™ d , K is related to II 
and V by 

k- 1 = v + it 1 (8) 

Thus, the matrix K describes the response of the system to an externally applied vector 
potential, including both the Hartree and non-Hartree parts of the interaction. In terms 
of a diagrammatic expansion, the function II is the fermion-hole proper polarization part 
(irreducible with respect to a single V line), while the function K is a sum of a geometric 
series, i.e., of a chain of proper polarization parts connected by single interaction lines V 
(see Fig. ©). 

The difference between the functions II and K should be particularly stressed in the 
context of the compressibility. Generally, the compressibility is given by the u = 0, q — > 
limit of the density linear response function. Denoting the ground state energy per unit 



volume by eo(n) where n is the density of particles, the inverse compressibility is ^ = n 2 ^^ 
W% - If the derivatives are taken with the external magnetic field B held fixed, the change 
in density induces a uniform magnetic field AB on the fermions. The fermions respond to 
the combination of the driving force and the induced uniform AB, and their response is 
therefore described by the linear response function K. The function II describes the case 
in which a change in the density is accompanied by a change in the external magnetic field 
B, in such a way that AB is kept constant. The two conditions yield two very different 
compressibilities, as can be seen by analyzing the case of non-interacting electrons around 
v — 1/2. If the density is changed around v = 1/2 with B being kept constant, - = 0. 



However, if AB is to be kept constant, the variation of the density results in a variation 
of the cyclotron frequency, and thus a variation of the kinetic energy of the lowest Landau 
level. Then, ! = We come back to this distinction in Sec. (5), where we discuss the 
gaps in the fractional quantum Hall states. 

Another use of Eq. (|6]) is for the analysis of the excitation modes of the system. The 
solutions of the equation in the absence of F are poles of n. They are characterized by a 
phase velocity u/q and by an eigenfunction v{6). Those solutions for which = J q>a; = 
are unaffected by the matrix V. They are poles of K as well. Solutions that do not satisfy 
these requirements are not poles of K. Some insight into the solutions of Eq. @ in the 
absence of a driving force is obtained by describing f(8) by its Fourier components, the 
"Landau parameters" 

fi = ^~ f T d6f(6)e ad (9) 

In "well behaved" Fermi liquids, the ratio m* /m is a finite number, and only a small, finite 
number of Landau parameters are appreciably different from zero. In that case Eq. (^) 
has a continuum of excitation modes at u < v F q, and a small number of discrete states at 
uj > v* F q. As we show in the next sections, the Fermi liquid we consider does not conform 
to this description. 

To conclude this section, we emphasize again the main points in the Fermi liquid picture 
we construct for the composite fermions. The energy of an excited state is described by 
the functional (0). The parameters in this functional, m* and /, will be calculated in the 
next section. The Landau function / does not include the Hartree contribution of either the 
Coulomb or the Chern-Simons interaction. 

The linear response of the system can be described in terms of two types of response 
functions. The first, n^(q, a>), describes the response of the function to the total driving 
force, while the second, if^q, u), describes the response of the system to the externally 
applied driving force. In the limit uo = and q — > 0, n describes the static compressibility 
when AB is kept fixed, while K describes the compressibility when B is kept fixed. The 
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poles of the linear response functions can be found, within Fermi liquid theory, by use of 
the equation of motion for n(k), Eq. (|6|). 



III. PERTURB ATIVE CALCULATION OF THE FERMI LIQUID PARAMETERS 

A. Starting Point 

In this section we calculate perturbatively the parameters characterizing the Fermi liquid 
state at AB = 0, namely m* and /. Landau's theory of the Fermi liquid expresses m* and / 
in terms of the self energy and the interaction operator of fermions at the Fermi level. The 
latter, in turn, are expressed in terms of the propagators for the fermions and the gauge 
field. We start this section by writing the propagators of the fermions and the gauge field. 
Then, we use these propagators to calculate the self energy, and use the self energy to show 
that the effective mass diverges. Then, we turn to discuss the Landau function, and pay a 
special attention to keeping internal consistencies between the approximations used. The 
Landau function we find turns out to have unique properties, which we discuss in detail. 

The fermion propagator in space-time representation is defined, as usual, as 



struction operators are expressed in the Heisenberg picture. The gauge field propagator is 
defined in a similar way. For the case of AB = 0, both the Hamiltonian (^) governing the 
time evolution of the system and its ground state are invariant to space and time trans- 
lations. Therefore, G is a function of r — r' and t — t' only, and it is diagonal in Fourier 
energy-momentum representation. The same applies also to the gauge field propagator, 
and, obviously, to the bare propagators. 

The bare propagator for the fermions is, for AB = 0, 




(10) 



where |0) is the ground state, t is the time ordering operator and the creation and de- 




1 



(11) 
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where \x is the chemical potential for the transformed fermions, and k is a wavevector. 

Since V ■ a = the gauge field propagator is a 2 x 2 matrix D^ u , with fi = denoting 
the time component and \i = 1 denoting the transverse spatial component. The gauge field 
propagator we use in the following calculation is obtained by the random phase approxima- 
tion (RPA). Its value for v = 1/0 (p = oo) was calculated in Ref. |1J, and was found to 
be 



/ V 1 

in „• g \ 



— -1-3- 

2vr 2tt. 



a 1 j *>L<±L n z ( 1 I i 

\ 2tt<P k F q "I Vi 27 rm ^ 27r0 2 eg'' / 



(12) 



In terms of the self energy of states at the Fermi level the effective mass m* is [14}] : 



j 8T, | 

1 _l_ o-^ 

^ d(, k \ k=k F 



m* = m ^, ° (13) 



where f A = |J - 



B. Lowest Order In — D approximation 

1. Fermion Effective Mass 

The self energy of the fermions is given, to first order in the gauge field propagator D, 
by the diagram given in Fig. 



^ = ^lS D <l k - k V - n W «) (M) 



(2vr 



There are two singularities in the evaluation of the integral flU). The first results from 
the high frequency regime of the longitudinal gauge propagator D 0Q . It does not affect the 
effective mass, and its discussion is therefore deferred to section (6). The second results from 
the infrared regime of _Dn, and is the source of the logarithmic divergence of the effective 
mass, as we now turn to discuss. In the long wavelength low frequency limit D is dominated 
by its diagonal transverse component, given by D u (q, cu) « D'f^ — g ^^J " 1 . We focus now 
on that element of D, and neglect all the others. The evaluation of the integrals in (13) is 
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similar to that carried out in studies of the electron-phonon coupling [15] . It is convenient to 
evaluate ui) — 0) and to introduce a variable q — |k — k'|. Most of the contribution 
to the k! integral comes from k' w k F . From the form of the gauge field propagator we 
observe that the important contribution to the q integral comes from q > 27r0^/|^^. The 
integral measure is transformed to 

roc r2ir roo rk+k' 2A"'fl 

/ k'dk' / d9 — >2 dk' dq . (15) 

J0 JO Jo J\k-k'\ J„2_(l c _V\2l(l,+ V\2_ n 2 



\k-k'\ y J q 2_( k _ k ty^( k + k jy_ q 

It is convenient to consider two regimes. In the first regime u > f k j F ^ k ~A F ^ , the integral 
measure can be approximated by 2 J °° dk' J 2kF dq. Then, the integral over k! restricts to 
satisfy < Q < u>, and the self energy becomes 

w , v 4> 2 ek F ( Ae 2 k F \ ,0 2 ek F 

L,(k,uj) — nOom — +2 -w (lb) 

V ' ; 2tt me 2 \ ^eu j A me 2 K ' 

Near the Fermi level the real part of E is larger than its imaginary part, thus justifying 
the notion of a quasi-particle. The origin of the logarithmic singularity is the ^ divergence 
of D for io = 0. Note that the important contributions to S come from a limited range 



of frequencies < f2 < to, but from a large range of wave- vectors 2^^^ < q < 2k F . 
Within the approximation used, £ is independent of k. Corrections to that approximation 
(e.g., to the approximation of the integral's measure) yield a weak and regular dependence 



on k. Substituting the approximation (|Tb|) in Eq. (|T3|) we get a logarithmically diverging 

duj ■ 



effective mass, Eq. (|2|). Note also that the residue (1 — |S) 1 vanishes at the Fermi level. 



In the second regime to < f 1 -^ — , the integral measure fll5|) is multiplied by a factor 
of k \ p and the singularity in the self energy is weaker. Near the Fermi surface the mass 
shell is defined by wlnu; cx \k — k F \. Thus, the first regime determines the quasi-particles 
dispersion relation and the effective mass is infinite. Later we shall argue that the above 
conclusions do not change when the set of diagrams included in the approximation for the 
self energy is extended. 



13 



2. Landau Interaction Function 



We now turn to discuss the Landau function /. For a rotationally invariant system the 
Landau function /(k, k') is a function of 9, the angle between the vectors k, k' (note that k « 
k' « fcjr). Standard Fermi liquid theory expresses the Landau function in terms of the Green's 
function residue z and the proper part of the interaction operator r(k 1; u\ \ k 2 , u 2 ] ki + q, ui + 
u; k 2 — q, u 2 — oo) as fI5|, 



/(k, k') = ^lim^of e (k, 0; k'0; k', u; k, -u) (17) 

where z^ F is the Green's function residue for poles on the Fermi surface, and 
r e (k, uj; k.'Q; k', Q; k, uj) = lim e ^ohm g ^or(k, u; k.'Q; k' + q, Q + e; k — q, uj — e). The full in- 
teraction operator T(k, k'Q; k' + q, Q + e; k — q, uj — e), including both the proper and 
improper parts, is defined through the two particle Green's function, in a way described in 
Fig. (^). The proper part, T, is the sum of all diagrams of T that are irreducible with respect 
to cutting a single interaction line, carrying a momentum q. The exclusion of the improper 
diagrams from the Landau function stems from the distinction between the linear response 
function II and K, defined in the previous section. The Landau function determines the 
linear response function II, through Eq. (|6p. Thus, it should include only diagrams that 
are "building blocks" for building the particle-hole proper polarization part (see Fig. ([I])). 
The diagrams excluded from the Landau function are those that serve to build the linear 
response function K. Examples are given in Fig. @). A comprehensive discussion of the 



distinction between V and V is given in Ref. |L7 



For our purpose, we find it convenient to rewrite Eq. (|T7D in a slightly different form. 
We define 

{G 2 (k, u)} € = lim^ G / (k, u + e)G/(k, uj) 

(18) 

{G 2 (k, cu)} q = lim^oG^k + q, u)G f (k, uj) 
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Since {G 2 (k,^)} e -{G 2 (k,w)} q = kp " F 6(u)6(k - k F ), ||- 0, Eq. (0) can be written 



as 

/(k,kO =lim_ ^j- f^(k,u;;k^;k'^ ; k,u;) [{G 2 (k, fi)} e - {G 2 (k, fi)} q 

(19) 



The extraction of an approximation for the Landau function from an approximation for 
the interaction operator and the Green's function residue should be done carefully, since the 
approximations for z^ F and for m* do not result from a systematic expansion in powers of a 
small parameter. Consistency between the approximations used is established by verifying 
that they satisfy the symmetries of the problem, namely, gauge invariance and galilean 
invariance. To be consistent with both symmetries, the approximations used should satisfy 
Ward's identities. Gauge invariance dictates the identity |p~5| . 



(Generally, Ward's identities relate S to T, and not to f . However, since here there is no 
direct contribution to the self energy, the contribution of the improper diagrams of T to the 
right hand side of Eq. (p0|) vanishes). When the self energy is evaluated by Eq. (0), Ward's 
identity ( p0|) suggests the approximation of the interaction operator by 

r (k, lu; k'fi; k', Q; k, u) = -^^^D(\k - k'|, u - SI) (21) 

and of the Green's function by Gf = G? . It can be shown that this approximation is 
consistent with the other Ward identities, too. 

Ward's identities are a test of consistency of the approximations for E and V. Attempting 
to use the consistent approximations for £ and V to formulate consistent approximations 
for the effective mass and the Landau function, we are in need for a similar consistency test. 
This test is given by the Fermi liquid identities, 



da 2tt i / 

— ■ — + JO j 

(22) 



dn m* 



J_ = J_ i A 

m m* 27r 
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It can be shown that if the approximations for T, Gf and £ used in Eq. ([19]) (where the 
approximation for E is used to determine m*) satisfy the Ward identity (|20"D, the identities 
(p2|) are satisfied, too. We show that for the identity involving fo in Appendix A. 

Consequently, in using Eq. ([19]) to consistently approximate the Landau function, the 
Green's function would be approximated by the interaction operator by Eq. (^), and 
the effective mass evaluated using the self energy (|16|) . Due to the divergence of m* at the 



Fermi energy, the limit uo — > of Eq. (|19|) should be taken carefully. Keeping to on the right 
hand side of (|i~9l) small but finite, we may define a frequency dependent Landau function, 
fuj(0)- Using Eq. (fl4|) for the self energy, Eq. ( ETC]) for the interaction operator and G^ for 
the Green function we find 

MO) = ^f-cos 2 -D(2k F sin-,u) (23) 
mm*{uj) 2 2 

For small lu, f u (9) is strongly peaked around 8 = 0, and the number of appreciably non-zero 
Landau parameters is very large. At the limit of zero frequency, 

lim^oA,(0) = ^S(9) (24) 
m 

and all Landau parameters are equal. The finite frequency m* (u) , f u (9) should be used 
when the equation of motion (|[) is analyzed. We present a partial analysis of this equation 
in the next section, and defer part of it to a future publication. 



C. Self Consistent Green's Function Approximation and Beyond 

The remainder of this section is devoted to an extension of the set of diagrams included 
in our approximation for £ and T. This extension suggests the conclusion that our ap- 
proximation for m*(uj), obtained from Eq. (|T4|), is exact in the limit uj — > 0, while our 
approximations for the Landau functions are not. 

The approximation we have used so far can be improved by evaluating the self energy 
self-consistent ly, i.e., by solving the equation, 
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The self energy obtained by Eq. ( p5|) is a sum of all "rainbow" diagrams. The singular 
parts of the effective mass and the Landau function are unaffected by this improvement of 
the approximation. To see that, follow the same route taken in evaluating Eq. (|14|). The 
integral over k! is proportional to the ratio 2fc ^™ = ^ . The residue is vanishingly small and 
the mass is infinitely large, but their product is finite, and is unaffected by the perturbation. 
Carrying out the rest of the integrals, we then find that the singular part of the self energy 



obtained from Eq. fl25|) is identical to that obtained from (0). Consequently, so is the 
effective mass. As for the Landau function, by using the relation between the bare mass, the 
effective mass and fi (Eq. (^)), it is easy to see that fi is unaffected, too. The same can 
be shown for f(6): when the self energy is approximated by the self consistent expression 
(f25|), Ward's identity requires the approximation of the interaction operator by a sum of 
ladder diagrams (see Fig. (^|)) and the Green function by the self consistent solution of Eq. 
(|H), Gf ~ [Q — e(k') — £(&', fi)] -1 . When these values of Y and Gf are substituted in Eq. 
(H)) i s found to be identical to that given in Eq. (|^). We carry out this calculation 
in Appendix B. 

A further improvement of the approximation ( fill) f° r the self energy might, in principle, 
be achieved by replacing the interaction vertices by dressed ones. As far as the effective mass 
is concerned, this replacement may lead to a modification of the logarithmic singularity in 
the limit uo — > 0, to a modification of the prefactor of the logarithm or to no modification of 
this singularity. We now argue that the last is true, i.e., that the effective mass (fjj is exact 
in the limit u — > 0. The essence of our argument is the observation, discussed below Eq. 
(fjjl), that for a state with an energy u close to the Fermi level, the contributions to the self 
energy come from small energy transfers < Q < u, but from a wide range of momentum 
transfers qo < q < 2k f, where go oc y/u. The logarithmic divergence of the effective mass is 
a result of the long wavelength contribution, of order of \fu. Thus, the interaction vertex 
of significance is that in which both the momentum q and the energy u transferred to the 

17 



gauge field approach zero, with the ratio u/q ~ -Juj — > 0. 

As a first step in our argument, consider separating the gauge field propagator D into long 
wavelengths (small q) and short wavelengths (large q) parts. The divergence of the effective 
mass results from the long wavelength part. Now, imagine repeating the calculation of 



the self energy leading to Eq. (|TJ|) or (|25|), but with the fermions' Green's functions, the 
gauge field propagator and the interaction vertices renormalized by the short wavelength 
gauge fluctuations. Moreover, imagine that this renormalization is done exactly, i.e., that 
Gf, D and the vertices are dressed by all possible lines of short wavelengths interactions 
(see Fig. (f|) for some low order contributions for the dressing of the vertex). What would 
be the effect of this renormalization on the logarithmic singularity in the self energy? Let 
us see how it affects each of the ingredients of the diagram in Fig. (Q). First, consider the 
fermion's Green's function. The interaction of the fermions' with the short wavelength gauge 
fluctuation induces a regular self energy in the fermions' Green function. We denote that self 
energy by S > (/c, uj). Second, consider the vertices. The self energy E > (/c,a;) determines not 
only the renormalization of the Green's function, but also that of the vertices. The latter is 
determined by a Ward identity: in the limit u/q — > 0, the interaction vertex is multiplied by 



1 + |T7J . Lastly, the gauge field propagator is not likely to be renormalized by the short 
wavelength gauge fluctuations: it is determined by the electron density and the electron's 
charge, none of which is subject to renormalization. 

Armed with these observations, we may now calculate the self energy due to the long 
wavelengths fluctuations, with the short wavelengths fluctuations taken into account exactly. 
The diagrams in Fig. (|2|) include two vertices, and each is multiplied by 1 + The integral 
over the magnitude of the internal momentum vector yields a factor of m*z = (1 + 
Thus, altogether, the short wavelengths fluctuations multiply the self energy flOp, which we 
denote by S°, by 1 + . The effective mass then becomes, 



m =m j as> ( 26 ) 

At sufficiently low frequencies, however, ^> 1, and m* is independent of S > , i.e., the 



diverging part of the effective mass is not affected by the short wavelength gauge fluctuations. 
We conclude, then, that the diagram leading to the divergence of the effective mass through a 
small energy-momentum exchange with the gauge field is not renormalized by large energy- 
momentum exchange. This conclusion suggests that our expression for the diverging part 
of the effective mass (0) is exact in the limit of u — > 0. 

Does our method take into account all diagrams contributing to the leading divergence 
in the self energy? Given a self energy diagram, in our method we single out the gauge field 
propagator with the smallest momentum, q min (denoted by the dashed line in Fig. (f|)), 
and assume that all the momenta carried by the other gauge field propagators (denoted by 
the dotted lines in Fig. (||)) are much larger. We thus effectively exclude the portions of 
phase space where two or more wave- vectors are close to q m in- This exclusion should not 
affect the coefficient of the divergent portion of if the self energy has contributions from 
a large range of momentum exchange. This is indeed the case for a Coulomb interaction, 
where the divergence is only logarithmic, and the momentum exchange goes up to 2kp. 
The separation of wave- vector scales leading to Eq. ( |2B| ) is therefore specific to the case of 
Coulomb interactions. 

The separation of length scales we use is in the spirit of a renormalization group analysis, 
which has been explicitly introduced by Nayak and Wilczek, and justified in the Coulomb 
case || ||. Nayak and Wilczek have also obtained a logarithmic singularity in the scaling 
of the frequency vs. \k — kp\ for the fermion propagator, but they have not discussed the 
value of the coefficient or consequences for Fermi liquid theory [fT2f| . 

The situation with the Landau interaction function / is considerably more complicated 
than that for the effective mass. Due to the identity (0), if m* diverges at the Fermi 
level we must have f\ = — . The other interaction parameters will be renormalized by the 
short wavelength gauge fluctuations, which are not treated correctly in our approximation. 
Thus, we are not able to make precise predictions for their behavior. In the limit where 
the bare mass m is small, so that the cyclotron frequency is large compared to the scale of 
electron-electron interaction, we expect the ground state and low-energy excitations to be 
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confined to the electronic lowest Landau level. Therefore, one might expect the energy cost 
of small deformations of the composite fermions' Fermi surface, Eq. (|5]), to be determined 
only by electron-electron interactions, and to be independent of the cyclotron energy, i.e., 
of m. In actuality, we see two exceptions to that rule. The first is a galilean boost of the 
Fermi surface (e.g., n(k) oc k x ). In that case the energy cost is determined by /i, which is, 
indeed, m-dependent. The second is n(k) = no, a uniform expansion of the Fermi surface. 
In that case the energy cost is determined by fo, which is related to the composite fermions 
compressibility. As discussed in section (2), the latter depends on the bare mass even when 
m —>■ 0. Consequently, we find that in the limit of m — > the first two Landau parameters, 
fo and fi, should have a value of the order of -jr. We believe, however, that for / > 2, 
the value of // is determined by the energy scale of electron-electron interaction. We also 
believe that lim/^oo/z ^ 0, so that f(Q) should include a delta function component, whose 
amplitude is determined by electron-electron interactions. At this stage, however, we are 
not able to substantiate this argument by a detailed calculation. 

In view of our caveats regarding the perturbative calculation of the Landau function, 
one may ask is there a parameter whose smallness could make our approximations for £ 
and T close to their exact values. This question is not easy to answer, since our results do 
not constitute the first few terms in a power expansion in a small parameter — The use of 
an RPA gauge field propagator involves a summation of all orders of perturbation theory. 
A small parameter in the sense defined above might be a parameter whose smallness makes 
the self energy ([16]) small. An obvious candidate is since lim^ S = 0. However, our 
results do not become exact when = 0. At = the magnetic field B = 2n(f)n = 0, 
and the problem is that of interacting electrons in the absence of a magnetic field. In the 
latter, which is of course a well studied problem, the electron-electron interaction gives rise 
to a wave-vector dependent self energy, which renormalizes the mass, in contrary to the 
vanishing self energy in our calculation for = 0. Consequently, our unperturbed problem 
could not be that of = and finite £ ^ s -, and has to be that of non-interacting electron 
at zero magnetic field, i.e., = e = 0. Our perturbative analysis should therefore apply 
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to the limit where both 0-^0 and — ► 0. For the self energy to be small, thus 
making the perturbative treatment sensible, the ratio between these two parameters should 
satisfy tjjz -> 0. In view of that, the approximation we have obtained for the Landau 
function can be better understood. At any finite non-integer 0, the problem we deal with 
is that of anyons in a magnetic field. Due to the anyonic statistics, the ground state of that 
problem is not confined to the lowest Landau level even as m — ► ||18|| . Thus, the bare mass 
affects the Landau function of the transformed fermions. It is only at the particular values 
of integer that when m — > the ground state is purely composed of lowest Landau level 
wave functions. Our perturbative treatment of the Landau function is not good enough to 
capture that feature. 

An alternative procedure which has been used to investigate a related problem of particles 
interacting with a gauge field is to introduce N species of fermions, and to analyze the 
problem in the limit of large N U ||. The RPA gauge field propagator in this model 
is identical to our Du, and the interaction with the gauge field induces a current-current 
interaction between the fermions. The leading term in the large N limit of that model are 
similar to the ones included in our self-consistent Green's function approximation. This 
large N procedure has not been applied directly to the Chern-Simons problem discussed 
here, however. 

In conclusion, our perturbative calculation indicates that the effective mass of the com- 
posite fermions diverges logarithmically, according to Eq. (Q), as the Fermi surface is ap- 
proached, while the Landau function approaches a delta function in 9. We believe that the 
former is an exact statement, while the latter is only approximate. In the next two sections 
we study possible consequences of these singularities on observable quantities. 
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IV. LINEAR RESPONSE FUNCTIONS AND COLLECTIVE MODES FOR THE 

COMPOSITE FERMIONS 



Having derived approximate values for the effective mass and the Landau function, we 
now turn to analyze the fermions' collective modes, and their linear response to a driving 
force. The collective modes we refer to here are the solutions of Eq. with F = 0, i.e., 
they are fictitious modes which would occur if the Hartree interaction were removed. As 
explained in section (2), some of these modes are unaffected by the inclusion of the Hartree 
term. Generally, in the absence of a driving force equation (|5]) has a continuum of solutions 
for which uj/q < v F , and might have a discrete spectrum of solutions for which uj/q > v* F . For 
any particular Landau function, the number of modes in the discrete part of the spectrum is 
smaller or equal to the number of non-zero Landau parameters. Within the framework used 
in the previous section we have found that the Landau function of the composite fermions is 
f{9 — 9') = ^-5(9 — 9'), i.e., its Landau parameters are all non-zero and equal. The effect 
of this particular Landau function turns out to be a trivial one: consider the slightly more 
general case of f(9 — 9') = 2nf5(9 — 9'). Substituting this Landau function in Eq. (|jD we find 
that its only effect is to shift the Fermi velocity from ^ to kp(^ + ^-). Galilean invariance, 
however, relates / and m* and requires ^* + 4^ = m IDJ- Thus, we find that the only effect 
of any delta function f{9) is to cancel the renormalization of the Fermi velocity and shift 
it back to its bare value. In the particular case we consider, m* = oo and f{9) = ^^-5(9). 
The Landau function shifts the mass back from infinity to its bare value. The effective mass 
m* then cancels from Eq. (f|), and, consequently, it does not affect the long wavelength low 
frequency behavior of linear response function. This result is of course consistent with the 
analyses of Kim et. al. || and Altshuler et. al. § for the response functions. 

The above results could be viewed as a limiting case of the behavior of the discrete 
spectrum when the number of appreciably non-zero Landau parameters is very large: the 
larger this number is, the larger is the number of modes at the range uj/q > v F , and the 
fuzzier is the distinction between the continuum and the discrete spectra. When all Landau 
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parameters are appreciably non-zero, the number of modes at v* F < u/q < vp becomes 
infinite, and the discrete spectrum may become continuous. 

Our use of the equation of motion (||) for a case in which the effective mass diverges at the 
Fermi level might be too crude an approximation. One may expect a more careful analysis 
to introduce a frequency dependence to m* and / in that equation, as presented in Eq. ( P5| ) . 
However, at low frequencies, the main feature we stressed here should not change. The 
Landau function has many appreciably non-zero Landau parameters, and thus the number 
of discrete collective modes is very large. At finite temperature or finite wavelengths these 
modes acquire a width, due to quasi-particle scattering, and the discrete spectrum becomes 
more of a continuum. 

We have seen that the singularities of m* and / mutually cancel in even denominator 
q, io — > linear response functions at zero temperature. This might cast some doubt on 
the physical significance of the divergence of the effective mass. Are there any observable 
quantities in which this divergence is manifested? The next section is devoted to the study 
of one such manifestation, in the variation of the chemical potential with the density at and 
between fractional quantum Hall states. 

V. THE VARIATION OF THE CHEMICAL POTENTIAL WITH THE DENSITY 
AT AND BETWEEN FRACTIONAL QUANTUM HALL STATES 

We now turn to consider the case of a large, yet finite, p, i.e., electrons at fractional 
quantized Hall states close to 1/0. The question to be discussed is the variation of the 
transformed fermions chemical potential \i as a function of the density n. For non-interacting 
electrons in a constant magnetic field, the chemical potential jumps discontinuously when 
the density is varied. The jumps take place at integer filling factors, and their height, Hu> c , 
is the energy gap for the integer quantized Hall effect states. For the problem we consider, 
care should be taken in defining an analogous question, since the density can be varied with 
the external magnetic field B being kept fixed, or with the effective magnetic field AB being 
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kept fixed. The different significance of the two quantities deserves further elaboration. 
The chemical potential p, is the energy cost involved in adding a transformed fermion to 
the system (along with a compensating contribution to the charge of the uniform positive 
background), while AB is held fixed. This chemical potential should be distinguished from 
the chemical potential p e for electrons, defined as the change in the ground state energy 
when one electron is added to the system and the external field B is held fixed. Since the 
total number of transformed fermions equals the total number of electrons iV e , we can simply 
write, 



,, — 9E I 
f 1 ~ dK\ AB 



.. _ dE I 

^ e dN e I B 



(27) 



where E is the system's ground state energy. 

Suppose that a fermion is added to the ground state at v — by use of the trans- 

formed fermion creation operator x I /+ (r), and suppose that the resulting wavefunction is 
projected onto the lowest available energy state. The added fermion carries a flux tube of 
<ft flux quanta, in a direction opposite to the external field B. If the latter is kept constant, 
this flux tube is uncompensated and the total flux becomes smaller by flux quanta. Thus, 
<fr fermions are pushed out of any one of the p occupied Landau levels, and the p + l'th 
level receives a total of 4>p + 1 fermions. If, however, it is AB that is kept constant, the 
flux tubes carried by the fermion are compensated by a change in B. Consequently, the 
total flux does not change, and only one fermion resides in the p + 1 Landau level. In that 
case, the resulting state has a fractional charge — e* in the vicinity of the point r, where 
e* = j^i, while the remaining charge — (e — e*) is pushed to the boundary of the system; 
i.e., a quasi-particle is added to the interior of the system. It follows from this analysis that 
the energy gap E g (v), defined as the energy to add one quasi-particle and one quasi-hole 
(well separated from each other) to the ground state at v — may be identified with a 
discontinuity in the fermion chemical potential \x at filling factor v. 

Our approach in studying /i(n) for a fixed AB is based on a perturbative calculation of 
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the energy cost involved in adding one fermion to the system in three different situations: 
(a) when p — 1 levels are completely filled, and the added fermion is the first one in the 
p'th level; (b) when the p'th level is almost filled, and the added fermion occupies the last 
vacant state on that level; (c) when p levels are completely filled and the added fermion is 
the first one in the p + l'th level. The three are denoted fi a {p) > A*&(p) > A*c(p) > respectively 
Obviously, // a (p) = t^cip ~ 1)- In the absence of the perturbation, for any given p, fi a = fib, 
and fi c — fib = Aui c = — . In the presence of the perturbation the jump in the chemical 
potential, \i c — fib, will be denoted by Au*. 

Our calculation of the three chemical potential differences is based on a perturbative 
study of the fermion propagator, defined in Eq. ([101) . One comment regarding the exact 
propagator is in place before we turn to perturbation theory. Since AB ^ 0, the Hamiltonian 
is not invariant to translations, and the fermion propagator is not diagonal in momentum 
representation. As we show below, however, it is diagonal in Landau level representation. 

The fermion propagator in Landau level representation is defined by 

G(n',m',t;n,m,t') = J ' d 2 r J m ,(r')G(r, t, r', f)0 n , m (r) (28) 

where (j) n ,m( r ) is the wavefunction for the state with angular momentum m and Landau level 
index n, and a symmetric gauge is used. As we now turn to show, G(n', m'; n, m) is non-zero 



only for m = m! and n — n'. The argument for that is based on a spherical geometry [19 
As discussed by Haldane, the problem of electrons in a magnetic field on a plane can be 
studied by using a spherical geometry and taking the limit of an infinite sphere. Similar 
to the plane, the single particle eigenstates of electrons on a sphere, subject to a magnetic 
field perpendicular to the sphere, are classified to energy-degenerate Landau levels. Unlike 
the plane, all states in the n'th Landau level are eigenstates of the angular momentum 
operator L 2 , with the eigenvalue (S + n)(S + n + 1), where 2S is the number of flux quanta 
threading the sphere. States within a Landau level can be taken to be eigenstates of L z , the 
^-component of the angular momentum, with integer eigenvalues ranging between — S — n 
to S + n. The spherical analog to a translationally invariant state is a state of zero total 
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angular momentum. As usual, G(n',m',t;n,m,t') is the amplitude of evolution from an 
initial state to a final state. Here, one state is a state in which a particle (or a hole) with 
angular momentum S + n and a ^-component m is added to a system that had previously 
no angular momentum. The other state is a state in which a particle (or a hole) with an 
angular momentum S + n' and a ^-component w! is added to a system that had previously 
no angular momentum. Both states are eigenstates of the angular momentum operators L 2 
and L z . Since the Hamiltonian conserves angular momentum, the amplitude for evolution 
from one state to the other is zero unless / = I 1 and m — m 1 . 

We shall begin with a calculation restricted to a lowest-order-in-D approximation, anal- 
ogous to that employed in section (3b) for the study of the effective mass m*. We shall see 
that in the limit of large p this calculation follows closely the calculations of the self energy 
and the effective mass for AB = 0. In particular, similar to the AB = case, the self en- 
ergy at the large p case is predominantly a function of frequency. At the end of the section 
we shall discuss corrections to the lowest-order-in-D approximation, and argue that the 
asymptotic form of the energy gap Au* for large p is unaffected by any corrections. 

The effect of the perturbation is studied through the self energy E. In the first-order- 
in-D approximation, E is calculated by use of the free fermion propagator, given by, 



where Aw c = AB/m, m is the bare mass, p is the chemical potential for the transformed 
fermions, and I is an index for the Landau level of the transformed fermion. 

The self energy of a state in the /'th Landau level is, in principle, a function of frequency, 
of the chemical potential (i.e., of the number of filled levels) and of I. It is given, in our 
first-order-in-/} approximation, by 



In the above equation, p is a chemical potential such that p levels are filled, tj = (j — l/2)Acu c 




(29) 




Q + uj + p — €j + ii]sgn(Q + uj) 



1 



D(q,Q) (30) 



and 
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M(l,j;q) = —J2\(l,g\J x (q)\j,g')\ 2 (31) 

N * 9,9' 

where is the Landau level degeneracy, j is a Landau level index, g, g' are indices of states 
within the Landau level and J± is the transverse current operator. 

In the following calculation of the self energy E we approximate the gauge field propagator 
D(q, uj) in ( p0[ ) by Eq. (|1^), i.e., by its value at p = oo. We denote the latter by D°°. We 
believe that this is a good approximation since for a large p most of the contribution to the q 
integral comes from q > R~ l , where R c is the cyclotron radius corresponding to the motion 
of a composite fermion in the p'th Landau level in a magnetic field AB. In that range of 
wave-vectors, the non-interacting fermions' response function (the fermions' free bubble) 
is relatively insensitive to the magnetic field. Since the dependence of the RPA gauge field 
propagator on p is only through its dependence on the fermions' response function, we expect 
the gauge field RPA propagator to be insensitive to AB, too. 

In terms of the self energy ( p0| ) . 

fi c - Hb = Au c + E(p + 1, Au* c , ii b {p)) - E(p, 0, fibip)) (32a) 
fJ>c- fJ>a = + E(p + 1, Au* c ,fi b (p)) - E(p, Au* c , fi b {p - 1)) (32b) 



fi b - fi a = E(p, 0, /Jtbip)) - E(p, Au*, /j, b (p - 1)) (32c) 
We start by calculating the energy gap Ao;* = \i c — /if,, and, in particular, its dependence 



on p. The two self energies appearing in Eq. ( |32a| ) differ in their Landau level index and 
frequency arguments, but have the same number of filled Landau levels. The self energy 
depends on the Landau level index only through the matrix elements, and this dependence 
is weak. We can therefore approximate 

A^ c * « Au, c (l + g) + Ef=i if f / $frM(p,r, q)Auj* c 

(33) 

x I I D°°(a Q) 
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In the limit p — > oo, approaches its value in the case of A73 = 0. As we saw in 
Section (3b), in our calculation of m* in the AB = case, this term could be neglected 
within the lowest-order-in-TJ approximation, since the a;-dependence of the self energy is 
more singular when uo — > 0. We make the same approximation here for large p. 

Since the gauge field propagator _D°°(q, Q) depends smoothly on Q and the matrix el- 
ements M.(p,j;q) depend smoothly on j, we can replace the sum over j by an integral 
■^j- J dej. Most of the contribution comes from €j ~ /i, and the limits of integration can be 
extended to include all the real axis. Carrying out this integration, we get, 

A ^ = ^ + i / If C f« M ^ ") < 34) 

where j = ~ p. 

Thus, the next step is to calculate M.(p,p;q). Assuming that q\\x, J±(q) = v y e iqx . The 
operators v y ,x can be expressed in terms of inter Landau level creation and destruction 
operators aj, a d and intra Landau levels creation and destruction operators 

in 



x = 75^(4 + a d + a+ + a g ) 



where Za = 1/vAB is the magnetic length corresponding to a magnetic field AB. In terms 
of these operators, 



M(p,p; q) 



^ q {W^ {a+d+ad) \p){ 



(36) 



The sum over g, g' cancels the in the denominator of Eq. (|3l|) 



Eq. (P6|) contains a derivative with respect to g of an expectation value of an operator 



at an harmonic oscillator energy eigenstate. Using R c = Ia\/P~ 1/2, we get, in the limit 
p -> oo, 

<p|e ig » = — ^ / 2 ^ dte^csinA^t = j Q ( gjRc ) (37) 

27T JO 

where Jo is Bessel's function. Consequently, 
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(Atu c R c ) 2 («-f) 2 for qR c « 1 



M(p,p; q) = {kuj c R c Ji{qR c )) 



2 



(AuJcRc) 2 ^^ cos 2 (qR c - 0.75tt) for 1 < qR c < R c /Ia 



negligible for g^A S> 1 

(38) 

Substituting into (U), and observing that in the large p limit most of the contribution comes 
from the intermediate wave-vector regime, we find 



-i 



A < dVL 1 



Auj* « Au c + v F / rfg / £>(q,fi) (39) 

JiJ" 1 JO 27T 7T 

where we used Au c R c = vf- Thus, the energy gap for the large p limit is, 

. */ \ e 2 kp7i 1 . . 

A ^ (P) = ^ P+1 )ln(2 P+ l) < 4 °» 

This expression for the energy gap is equivalent to the one given in the introduction, Eq. 
(|) (note that for large p, In (2p + 1) ~ In (4>p + 1) ~ hip, and k F ~ /^ 1 (2/0) 1 / 2 ). It is 
also consistent with the recent calculation of Kim et. al. [^] . The origin of the logarithmic 



singularity is the same as its origin in the calculation of the effective mass for the v = l/<p 
case, namely, the ^ divergence of D°°(q, 0). In fact, after making the approximation of the 
sum over j by an integral (Eq. ([34])), we brought the self energy in the large p case to a 
form very similar to that of the self energy in the AB = case, with the main difference 
being the cut-off of the infra-red divergence in the momentum integral, resulting from the 
suppression of the matrix elements for q -C R~ l . Thus, the magnitude of the discontinuous 
jump in the chemical potential at the fractional quantized Hall states reflects the logarithmic 
singularity of the fermions' effective mass. For any finite p, this singularity is cut-off by 
lnf^ ~ ln(2p + 1). 

Having calculated the energy gap, we now turn to fi c {p) — yUa(p)- The sum of matrix 
elements, A4(l,j;q), depends strongly on \l — j\, and only weakly on l,j independently. 
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Consequently, the self energy difference appearing in Eq. ( |32b|) is regular, and fi c {p)—fi a {p) ~ 
Au c . Similarly, fi b (p) - fi b (p - 1) « Au c . 

Lastly, we turn to //& — p, a , the last of the three chemical potential differences defined in 



Eqs. fl52D- Obviously, it is uniquely determined by Eqs. ( |32aj ) and (|32b[ ). It is, however, 
interesting to carry out an independent calculation of this difference in chemical potentials. 
Unlike the self energy difference \i c — fi b , the self energy difference appearing in Eq. ( |32c|) 
involves self energies at different chemical potentials, i.e., different number of filled levels. 
It is given by, 



(2tt)2 



Sl+/j, b {p)—ej+irj Cl+Ab>*+Hb(p— 1)— ej+iri 



- Ef= P+ i M (p, j ; g) J f ( ^tw^n+A.g^^-^-e.-iJ ^Cq^ ^) (41) 



§ (/it - /ia) + / t|&A1(p,P; g)£>(q, e p - A< - // 6 (p - I)) 



(2^)2 

In writing this expression we have separated the difference between the two self energies to 
a part that does not involve a change in the level occupation and a part that describes the 
change of the p'th level from vacant to occupied. The first part varies smoothly when the 
frequency and chemical potential are varied, and is therefore approximated by a derivative. 
In the second part, in the large p limit, D(q,e p — Auj* — fi b (p — l)) can be replaced by D(q, 0). 
Following these approximations, we indeed find that 

^-^f ln(2p+l) 
f^ln(2p+l) 

(42) 



/i6 - /ia = AUJ C 171 



= Auj c - Alu* 

Eqs. fl40D and (^) confirm the picture of /i(n) for fixed AB given in Fig. (||): the 
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chemical potential jumps discontinuously whenever the fermions start occupying a previously 
vacant level. The magnitude of that jump is singularly reduced by the interaction with the 
gauge field fluctuations. In between two discontinuous jumps, the chemical potential is not 
constant. The total change of the chemical potential during the population of a level is 
Aui c — Au*, and it approaches Auj c as p — > oo, i.e., as the discontinuous jump is suppressed. 
Our calculation does not give any information about the way the chemical potential varies 
between p, a and One could imagine a continuous linear variation, which would result 
from a Hartree-Fock approximation, or a discontinuous variation, due to, e.g., the formation 
of fractional quantum Hall states of the composite fermions. Considering the fact that the 
fractional quantum Hall effect has only been observed in the first few Landau levels, and 
the fact that states in high Landau levels are extended over a radius of R c ^> we believe 
that a linear variation is at least an excellent approximation. 

Thus far, we have restricted our analysis to a lowest-order-in-D approximation for the 
self energy at AB ^ 0. We found that the energy gap Au* obtained in this approximation 
is consistent with Eq. (JJ) and with the logarithmically diverging mass obtained by a similar 
approximation for the AB = case (see Section (3b)). In Section (3c) we argued that the 
calculated logarithmic divergence of m*, together with its coefficient, are in fact exact in the 
oo — > limit. This is because a Ward identity guarantees that corrections to the vertices due 
to short wavelength fluctuations are cancelled by the correction to the intermediate fermion 
Green's function (see Fig. (f|)) and by the factor (1 + which enters the expression 

for m* (Eq. (|13"D). A similar cancelation takes place when we extend our calculation of the 
energy gap Au* beyond the lowest-order-in-D calculation. Again, we find that corrections 
to the vertices are cancelled by corrections to the intermediate fermion Green's function and 
by the factor (1 + |^-) neglected in Eq. (|3^). Thus, we find that Eq. ( ^0|) remains exact in 
the limit p — ► oo. Due to the similarity with the analysis presented in Section (3c), we omit 
the details here. 

To conclude this section, we make one more comment regarding the energy gap. Our 
calculation of the FQHE energy gaps is based on the study of poles of a single fermion 
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Green's function. An alternative way to extract the energy gap /S.u* is through the linear 
response functions II and K at an integer fermion's filling factor p. As illustrated by an 
RPA calculation |]J and a modified RPA calculation ||22|| , in the large wave-vector limit, 
q ^> 2kp, the lowest poles of both II and K occur at a frequency that equals the energy gap. 
These poles describe the excitation of a well separated quasi-hole quasi-particle pair. The 
modified RPA calculation employs an approximation in which only one Landau parameter, 
/i, is non-zero, and shows that while at the large wave- vector limit the excitation energies 
are unaffected by the Landau parameter, this is not the case for small values of q. In 
the latter case the quasi-particle and quasi-hole are not well separated, and their mutual 
interaction affects the energy of the excitation. 

VI. FERMION GREEN'S FUNCTION AND COUPLING TO D 00 

As has been noted above, the logarithmic divergence of |j, arising from interactions 
with the gauge field, implies that the weight of the quasi-particle contribution to the 
fermion Green's function Gf is proportional to l/ln|fc — kp\ when k — > kp, and vanishes 
at the Fermi surface. However, as originally noted by HLR, there is another effect which 
causes Gf to vanish for any value of u, for any momentum and any finite energy, in the limit 
of infinite system size. This effect arises from of the coupling of the fermion propagator to 
the longitudinal gauge propagator D Q0 , which thus far has been neglected in our discussion. 

The physical reason for this divergence is that the transformed fermion operator \l/ + (r) 
has the effect of adding an electron at point r, and instantaneously turning on a solenoid 
containing two flux quanta. The impulse electric field generated by the solenoid produces 
inter Landau level excitations at large distances from r (i.e., it produces long-wavelength 
magnetoplasma excitations at the cyclotron frequency u c ) whose mean number diverges 
logarithmically with the size of the system. The probability that the many electron state 
has remained entirely in the lowest Landau level (or in any other state with only a finite 
number of cyclotron excitations) is therefore found to vanish as a power of the size of the 



32 



system 0. 

In order to have a non vanishing weight for the one fermion Green's function at low en- 
ergies, it is necessary to define a renormalized version of the Green function. One possibility 
would be to carry out calculations for a large but finite system, and to renormalize the over- 
all weight by a size-dependent factor before taking the limit of infinite size. An alternative 
approach is to work in an infinite system in the limit where all operators are projected onto 
the lowest Landau level of the electrons. This is equivalent to taking the limit where the 
bare mass m — > (i.e., u c — > oo), while the electron-electron interaction is held constant. 
Simultaneously, the flux in the solenoid associated with \l/ + (r) should be turned on at a rate 
slow compared with uj c but fast compared to the scale of the electron-electron interaction. 

Although neither of these procedures has been carried out in detail, we expect that the 
renormalized fermion Green's function should indeed have the low energy properties obtained 
for Gf in the previous sections, where the coupling to Dqq was neglected. In particular, the 
renormalized Green function should have a singularity at u — (k — hp), whose weight 
vanishes as 1/ln \k — kp\. There is no indication that the coupling to Dqo would change 
the singular behavior of the quasi-particle dispersion for k — > kp, or of the energy gaps at 

f = T^"T- 

(f>P+l 

VII. SUMMARY 

In conclusion, in this work we discussed the effect of the diverging long wavelength low 
frequency fluctuations of the Chern— Simons gauge field on the properties of the electronic 
state at and near v — 1/0. 

In Sec. (2) we constructed an extension of Fermi liquid theory designed to describe 
the effect of the Chern-Simons fluctuations on the low energy excitations of the composite 
fermions. We emphasized the separation between the Hartree and non-Hartree part of the 
interaction. 

In Sec (3) we derived approximate expressions for the Fermi liquid parameters, namely 
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the fermions' effective mass and the Landau function. We showed that at v = 1/0 both 
behave singularly. The effective mass diverges as the Fermi surface is approached, and the 
Landau function has a part that approaches a delta function at 6 — 0. We argued that our 
approximation for the effective mass becomes exact in the limit u — > 0. 

In section (4) We showed that at AB = the singularities in m* and f(9) are cancelled 
in the linear response functions in the limit of low frequencies and long wavelengths. 

In sec. (5) we investigated the variation of the chemical potential of the transformed 
fermions when the density is varied and AB is kept fixed. We showed that within our 
perturbative treatment, the energy gap of the fractional Hall states depends singularly on 
the filling factor, but the overall slope of n(n) is not subject to that singularity. 

Finally, in Sec. (6) we briefly reviewed the effect of the fermions coupling to longitudinal 
gauge field fluctuations on their Green's function. We argued that this coupling does not 
affect the results obtained in the previous sections. 

Even within the limited context of clean systems, to which we have constrained ourselves 
in this paper, we left several important questions opened. We would like to point out 
some of these. First, while we believe our approximation for the effective mass of the 
fermions becomes exact at a certain limit, we do not have a similar understanding of our 
approximation for the Landau function. We consider it likely that our approximation for 
the Landau function yields the correct type of singularity (a delta function), but with an 
incorrect coefficient. 

Second, the analysis of this paper has been explicitly confined to the case of Coulomb 
interactions between the electrons (i.e., an interaction that behaves as 1/r at large separa- 
tions. Extensions to a shorter range interaction would be of considerable theoretical interest, 
but must be carried out with care. Perhaps most importantly, care should be exercised in 
discussing the very notions of an effective mass and Fermi liquid parameters, since the de- 
cay rate of the quasiparticle might be comparable to its energy (in contrast to the Coulomb 
interaction, where the real part of the energy is larger than the imaginary part by a factor 
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of In uo , a factor that diverges as the Fermi surface is approached) . 

Finally, it should also be noted that in most practical applications to present experiments 
on electron systems in partially filled Landau level, the logarithmic contribution to the 
FQHE energy gap is smaller than the finite, non-singular, contributions arising from short 
wavelength fluctuations in the gauge field, which we were not able to calculate reliably. 
These contributions depend on the behavior of the electron-electron interaction at short 
distances, and are therefore "non-universal". 
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APPENDIX A: CONSISTENCY IN THE APPROXIMATIONS FOR THE 
EFFECTIVE MASS AND THE LANDAU FUNCTION 

In this appendix we show that the approximation scheme we have used for approximat- 
ing the Landau function does satisfy the first of the Fermi liquid identities given in (f22|). 
Consistency with the second Ward identity can be shown using similar methods. The first 
identity is 

d f--, + h (Al) 
an m 



As described in section (3), we write the zeroth moment of the Landau function as [T3 
Q 

35 



/ = lim^o-^- f-^f ^f e (k,w;k'fi;k',fi;k,a;) [{G 2 (k, Q)} t - {G 2 (k, fi)} t 



(A2) 



where 



{G 2 (k, cu)} e = lim^ G/(k, w + e)G/(k, u) 

(A3) 

{G 2 (k,^)} q = lim q ^ G/(k + q s w)G / (k,u;) 

We would like to substitute approximate values for r e , G/, m* in Eq. (|A2|) , in order to obtain 
an approximation for fo. The choice of this approximate values, however, should be such 
that the identity ( |Al| ) is satisfied. We now turn to explain how this choice should be made. 

Our search for an approximation consistent with ( |Al| ) starts with expressing / in terms 
of derivatives of the self energy S. In doing that, we use a set of Ward identities, one of 



which was already given above We refer the reader to Ref. [151 for a comprehensive 

discussion of this set of identities. A shortened notation is useful here. We write Eq. (|A2|) as 
f = =2^ij[{G 2 } e -{G 2 } q ]f e . Ward's identity ©, which is useful here, is £ = ijf e {G 2 } e . 
Also useful are the identities Jl6|] [15[ 



r = f q -i / f e [{G 2 } e - {G 2 } q ]f q (A4) 



(with r 9 (k, lu; k'f2; k', Q; k, u) = lim^oh m e->or(k, tv; k'f2; k' + q, Q + e; k — q, uj — e)), and 
the Ward identity for the chemical potential, 

(A5, 

Substituting these three identities into Eq. ( |A2|) we get, 



leading to, 



(A6) 

_ 2?r 8S I 2tt <9S I 3Sf 
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Note that all derivatives should be taken at the Fermi level. 

Having expressed /o in terms of derivatives of the self energy, we now express j£ in similar 
terms. The self energy is, generally, a function of k and u (we deal with the AB = case), 
and it depends also on the chemical potential /x. To stress that, we write it as E(k, a>;//). 
The chemical potential is the energy needed for adding one fermion at the Fermi level. It is 
given by, 

\i = e kF + E(k F , 0; //) (A8) 

Taking the complete derivative |- of both sides of (|A8|) , we get, 

i i as i 

d/l _ 27T ~*~ de k \k = k F 



e fc 

dn m 1— 

Oil 



2, 1-iS 



(A9) 



771* 1— ^ 
OjJ, 



Comparing with Eq. ( [AT] ) we see that the approximate value for the effective mass 
used in Eq. (|A2|) to determine the Landau function should be the same as the one used 
in the identity (|A1|) . Moreover, we get to the following prescription for establishing an 



approximation for m* and / that respects the Fermi liquid identities (p2 ): 



The first step is to choose an approximation for the Green's function Gt , the interaction 
operator T and the self energy S in such a way that the Ward identities (20) and (A5) 
are satisfied. 

At the second step, the approximate value for S should be used to extract an approx- 
imation for the effective mass. 

Finally, the approximate values for the effective mass, the interaction operator and the 
Green's function should be substituted in Eq. (0), to obtain an approximate value 
for the Landau function. 
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The application of this prescription to the first-order-in-D approximation we have used 
in Section (3) is straightforward. As we discussed in Section (3) (below Eq. (p0|)) the Ward 
identity fl20|) is satisfied if the fermion Green function is approximated by the free one G®, 
the interaction operator by Eq. ([H]) and the self energy by (|i4|). Thus, in consistently 



approximating the Landau function using Eq. ( |A2[ ), these are the values that should be 
used. In particular, the effective mass used should be the logarithmically diverging one @. 
Following these principles, the approximation obtained for the Landau function is the zero 
frequency limit of Eq. (^). The application of the above prescription to the self-consistent 
Green's function approximation is slightly more complicated, and is discussed in Appendix 
B. 

To conclude, we note that the observations made above regarding the consistency of 
approximations for the Landau functions are hard to make using the conventional expression 
for the Landau function, Eq. ([[7]). The alternative equivalent expression we suggested, Eq. 
(|I9|), makes this analysis simpler. 

APPENDIX B: APPROXIMATION OF THE INTERACTION OPERATOR AS A 

SUM OF LADDER DIAGRAMS. 

In this appendix we show that the self consistent Green's function approximation for 



the self energy leads to the same Landau function, Eq. (|23D , as the one obtained from the 
first-order-in-D approximation. In the self consistent approximation, the Green's function 
is G « [fi - e(Jfe') - E(Jfe', VI)}- 1 = z(n)(Q - v* F (k - k F )Y\ where z(Q) is the residue. As we 
saw in section (3), the leading singularity of the self energy, given by (0), is the same as 
the one obtained in the first-order-in-D approximation. The interaction operator should 



be approximated in such a way that the Ward identity flzOD is satisfied. We start by finding 
out what this way is. 

Taking the derivative of (|25| ) with respect to frequency we get, 
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OS _ r d 2 k' r dtt 1 (kxk') 2 
dui J (2tt) 2 J 2tt m 2 jk-k'j 2 

(Bl) 

xD n (\k - k'|,Q) ( „_ n _ e(y) i E(fc , ih ,_ n) ) 2 (1 - §Un) 

Substituting recursively |j into the right hand side of ( pi]) and comparing with Eq. fl20|), we 

find that for Ward's identity to be satisfied, the interaction operator should be approximated 

by the sum of ladder diagrams, illustrated in Fig (|). 

Ward's identity requires |^ = i JT e {G 2 } t . Since the left hand side of the identity is 

unchanged when the first-order-in-/} approximation is changed to a self consistent one, so 

should be the case with the right hand side. In the self consistent approximation the two 

Green functions carry two powers of the residue. The integral over the magnitude of the 

intermediate momentum yields a factor of -V ~ rn* ~ z^ 1 . Therefore, for the right hand 

Up- 
side of the Ward identity to be independent of z, the interaction operator should be 0(^ _1 ). 

We now show that this observation is consistent with the equation defining the ladder sum 

for the interaction operator, and then turn to examine its implications on the approximation 

for the Landau function. 

The ladder diagram approximation for the interaction operator is the solution to the 

equation, 

r (k, to; k'fi; k, fi; k, to) = -^%^D(\k - k'|, u - Q) 



m 2 |k-k'| 2 

[ 

(2T) 2 J "2tT m 2 " |k-k"| 2 



+ / / ^^CTF^d k " k "U " uj")G 2 (k", uf)Y e (k", J'- kV; k', J- k", u") 



(B2) 



If the interaction operator is to be 0(z x ), the first term on the right hand side is negligible 
relative to the second. Neglecting that term, equation ( |B2| ) becomes an eigenvalue equation 



for f . For this equation to be consistent with our z-power counting, the second term of the 
right hand side should be of the same order of z as T. This is indeed the case: again, the 
two Green's functions carry two powers of the residue, and the integral over k" introduces 
a factor of (vp)^ 1 ~ m* ~ z~ l . Now, our interest is in momenta k, k' close to kp and 
frequencies u,Q close to zero. The denominators in D and G 2 constrain k" to be close to 

39 



kp and lu" 0. We may therefore approximate \k — k"\ — 2kp sin where 9" is the angle 
between k and k". The integral over 9" then diverges at small angles, due to the divergence 
of D. This logarithmic divergence is the same as the one leading to the divergence of m* 
and z~ x . Altogether, then, the right hand side of ( |B2| ) has the same order of z as the left 
one, namely, the order of T. We therefore conclude that T ~ 0(z~ l ). 

Finally, we use the same type of z-power counting to examine the approximation for 
the Landau function. We start with Eq. (Jl9|). Within the self-consistent Green's function 
approximation, the Green's functions substituted in this equation carry two powers of the 
residue z, and the interaction operator, as shown above, is of order z~ l . Again, the integral 
over k' introduces a factor of 4- ~ z~ x , and the Landau function becomes 0(z°), i.e., it 

V F 

is the same as the one obtained in the first-order-in-D approximation, the zero frequency 
limit of Eq. (||). 
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FIGURES 

k+q, e+co 
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FIG. 1. The diagrammatic representation of the linear response functions II and K. The 
response function II is the fully dressed fermion-hole bubble, irreducible with respect to a single 
gauge field line. The response function K is the geometric series summing strings of bubbles 
connected by single gauge field propagator line. Solid lines are fermion propagators, while dashed 
lines are gauge field propagators. The physical meaning of both is discussed in the text. 
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FIG. 2. 

The self energy diagrams, (a.) The lowest order in D diagram, leading to Eq. (0). The 
dashed line represents the RPA expression for the gauge field propagator (Eq. (|T2|) . (b.) 
The self consistent approximation, leading to Eq. fl25|). 
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FIG. 3. 

The proper interaction operator T and the approximations we employ for it. (a.) The 
interaction operator T is defined via the fermions' two particle Green function G^ 2 \ In 
this figure, solid lines represent exact single particle Green's functions, (b.) The simplest 
approximation for T, used in calculating Eq. (^l|). The dashed lines represent the RPA 
expression for the gauge field propagator (Eq. fll~2|)), and the crosses represent bare vertices. 
The left diagram is proper, and is included in T. The right diagram is improper, and is not 
included in T. (c.) The ladder diagrams sum for T, summed in Appendix B. 
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FIG. 4. 

(a.) Some low order contributions to the dressed vertex. The dashed line represents a 
low frequency long wavelength gauge field propagator. The dotted lines are short wavelength 
gauge field propagators, (b.) The renormalization of the self energy by short wavelength 
gauge field fluctuations. The dashed line carries the lowest momentum in the diagram. The 
double line (solid and dotted) denotes a fermion's Green's function, renormalized by the 
short wavelength gauge field fluctuations. As explained in text, the dominant contribution 
from the gauge field propagator is not expected to be renormalized. 
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FIG. 5. The variation of the chemical potential when the density is varied, and AB is held 
fixed. The location of the origin and the units are arbitrary. Discontinuous jumps take place 
when a vacant Landau level starts being populated. Their magnitudes are the energy gaps for 
the corresponding fractional quantum Hall states. They reflect the value of m* . Smooth variation 
takes place in between discontinuous jumps. The average slope of the curve is independent of the 
singular part of m* . 
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